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Abstract 

in 

' The PAMELA apparatus has been assembled and it is ready to be launched in a satellite 

^ ■ mission to study mainly the antiparticle component of cosmic rays. In this paper the per- 

I \ formances obtained for the silicon microstrip detectors used in the magnetic spectrometer 

I are presented. This subdetector reconstructs the curvature of a charged particle in the mag- 

• netic field produced by a permanent magnet and consequently determines momentum and 

charge sign, thanks to a very good accuracy in the position measurements (better than 3 ^ra 
\ in the bending coordinate). A complete simulation of the silicon microstrip detectors has 

' been developed in order to investigate in great detail the sensor's characteristics. Simu- 

^ lated events have been then compared with data gathered from minimum ionizing particle 

(MIP) beams during the last years in order to tune free parameters of the simulation. Finally 
some either widely used or original position finding algorithms, designed for such kind of 
detectors, have been applied to events with different incidence angles. As a result of the 
analysis, a method of impact point reconstruction can be chosen, depending on both the 
particle's incidence angle and the cluster multiplicity, so as to maximize the capability of 
the spectrometer in antiparticle tagging. 
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Fig. 1 . Schematic lateral view of the PAMELA detector (left) and a photograph of it 
(right) taken before the delivery of the instrument to the Russian space company. The 
geometrical acceptance of the detector is 20.5 cm^sr. 



1 Introduction 

The PAMELA telescope [1,2] will be put in orbit within the 2005 on board of the 
Resurs DKl Russian satellite for a three-year-long mission on a quasi-polar orbit 
(70.4 deg. inclination, 350 to 600 km height) to study the cosmic ray flux, with 
a special interest on the antimatter component. The detector is composed of sev- 
eral subsystems, schematically shown in fig. 1: a Time of Flight (ToF) apparatus, 
which also provides the trigger signal, a solid state Magnetic Spectrometer [3,4], 
surrounded by an anticoincidence shield, and an Electromagnetic Calorimeter [5] 
in which single-sided silicon detector planes are interleaved with tungsten absorber 
up to a total thickness of about 16 radiation lengths. Anticoincidence scintillators 
define the external geometry of the detector and their signals will be exploited in 
the off-line rejection of spurious tracks; below the calorimeter another scintillator 
plane (S4) and a Neutron Detector can provide additional information when show- 
ers are not fully contained in the calorimeter. 

The whole detector can discriminate protons and electrons from their antiparticles 
and provide energy measurements; also Z > 1 nuclei may cross the entire spec- 
trometer and consequently can be detected. Antimatter acquisitions will cover the 
energy range 80 MeV-190 GeV for antiprotons and 50 MeV-270 GeV for positrons, 
so extending the current upper limit (about 50 GeV) given by balloon-borne mis- 
sions [6]. In addition, the long duration of the flight will provide larger statistics 
in particular in the high-energy range, where the power-law spectrum of cosmic 
rays requires either large acceptances or long exposure times ^ . The obtained high- 

^ The differential flux of galactic cosmic rays at 1 Astronomical Unit from the Sun is 
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energy antiparticle flux wiU constrain the models describing the origin of antimatter 
in the Universe. The apparatus can be exploited also to find dark-matter signatures 
through the detection of high-energy antiprotons originating from neutralino anni- 
hilations [8]. 



2 The magnetic spectrometer 

2.1 Magnet 

The magnetic spectrometer is the core of the PAMELA apparatus: it is based on 
a permanent magnet and consists of six detector planes which measure both the 
impact coordinates of the incoming particles. The magnet is made of a Nd-Fe- 
B alloy, with a high value of the residual magnetization (about 1.3 T). Blocks of 
magnetic material define a rectangular cavity (132 x 162 mm^) where particles are 
detected. Here the magnetic field is roughly uniform and oriented along the Y co- 
ordinate in the PAMELA reference frame (fig. 1). As a consequence, particles are 
bent in the XZ plane within the cavity, due to the Lorentz force F = gv x B. Five 
identical magnetic modules, each one 80 mm high, are superimposed each other 
and interleaved with six detector planes, which are inserted in the magnetic tower 
by means of dedicated slits. The magnetic field in the centre of the cavity is 0.48 T. 
Measurements of the three components of the magnetic field have been taken at a 
fixed pitch in about 70 000 points inside the cavity. Such values will be used dur- 
ing the off-line analysis to precisely determine the particle's momentum through a 
numerical integration of its equation of motion inside the cavity. 

2.2 Silicon detectors 

When the characteristics of the PAMELA experiment were studied, the main re- 
quirements of the detectors to be used in the magnetic spectrometer were defined. 
Essentially they can be listed as in the following: 

a) provide two coordinates per detector; 

b) maximize the spatial resolution, in particular for the bending coordinate; 

c) minimize the multiple scattering. 

The best candidates to comply with such requirements appeared to be double-sided 
silicon microstrip detectors with a small strip pitch, associated with low-noise elec- 
tronics. The multiple scattering was limited by excluding any dead layer above or 
beneath the detectors. In the resulting configuration six detectors, manufactured by 

proportional to in the energy range ~ 10 — 10^ GeV/n [7]. 
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Fig. 2. Photograph of a siUcon detector plane. 

Hamamatsu Photonics [9], are arranged in each plane (fig. 2). A detector is made of 
a 300 /im thick n-type silicon wafer, with p-type strips implanted at 25.5 fim pitch 
on the junction side. These strips are used to determine the X coordinate of the 
incoming particle ^ . n"''-type strips, orthogonal to the previous ones, are implanted 
with 66.5 fim pitch on the Y (ohmic) side. The read-out pitch on the junction side 
is actually 51 fim, since a p strip out of two is left floating: charge collected by that 
strip influences the output of the adjacent read-out channels [10] . On the ohmic side 
a p blocking strip is added between two consecutive n strips in order to increase 
the interstrip resistance. Integrated decoupling capacitors are present on both sides: 
they are obtained by depositing an insulating layer (0.1 /im Si02) and metallic strips 
above the implants. On the ohmic view a second metal layer, separated by another 
Si02 deposition, 5 fim thick, brings the read-out metallic strips parallel to the junc- 
tion side's ones. Both surfaces of the detectors are passivated; electrical signals can 
be taken out from metallic pads located at both ends of the strips. The depletion 
voltage given by the manufacturer is 60 ± 20 V and 80 V bias will be supplied 
during the operational phase. The detectors are arranged in three modules named 
ladders: two sensors and a double-sided hybrid circuit form each ladder. These 
elements are glued in such a way that X side strips meet at their ends: electrical 
connections among the corresponding strips are provided by 17/Lim wire bonds on 
both sides. On the ohmic view two distinct strips, about 70 mm apart, are connected 
to the same electronic channel, so giving rise to a position ambiguity, that can be 
easily resolved on the basis of the other PAMELA'S subdetectors. The percentage 
of bad strips, due mainly to unconnected or noisy channels, ranges between 1% and 
5% for all the ladders. The mechanical structure of the plane is shown in fig. 2 and 
it minimizes the multiple scattering since no additional dead layer is present on the 
particle's trajectory. Nevertheless, in order to increase the rigidity of the structure. 



^ The X coordinate is related to the particle's curvature: for this reason the X view of the 
detector is also called the bending view. 
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Table 1 

Geometrical and electrical characteristics of the silicon microstrip detectors. 



four carbon-fibre bars (visible in the photograph) with very high Young modulus 
(about 300 GPa) are glued at both sides of each ladder. The resulting plane has the 
first resonance frequency located beyond 300 Hz and it has survived vibrational 
tests whose intensities were twice those expected during the launch of the satel- 
lite [11]. A summary of geometrical and electrical characteristics of the detectors 
is given in table 1 . 

2.3 Front-end electronics 

The front-end electronics is housed on the double-sided hybrid: it is based upon 
the VAl [12] chip, which consists of 128 charge preamplifiers, shapers and sample- 
and-hold stages in series with a multiplexer controlled by a shift register. The dy- 
namic range of the chip corresponds approximately to ±10 MIP's. The analog sig- 
nal on each channel is then digitized by means of a 12-bit ADC per ladder view. 
As known in literature [13], the noise observed in a detector-preamplifier chain is 
essentially related to four parameters, depending on both the sensitive element and 
the associated front-end electronics: the bias resistance Rb, the leakage current i^, 
the total capacitance C to ground seen by the preamplifier and the shaping time 
r. For a given shaping time (r = 1 //s in our case), the noise is minimized by a 
large value of the bias resistance and by small values of leakage current and load 
capacitance. Let us briefly discuss each contribution as it appears in our configu- 
ration. Bias resistors are used to connect each strip to the guard rings, which are 
kept at a fixed potential. On the junction side the punch-through technology [14] 
with an additional foxfet electrode [15] has been exploited to have Rb > 50 Ml]; 
on the ohmic side polysilicon resistors have been used and the corresponding re- 
sistance is greater than lOMfi. The measured leakage current is of the order of 
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1 nA/strip. The Rb value and the limit on iL have been designed to give a negligible 
noise, compared to the contribution from the capacitive load seen by the preampli- 
fier. The requirement on this parameter imposed to the detector's manufacturer was 
C < 10 pF (20 pF) for the junction (ohmic) side of a sensor. 

The VAl Equivalent Noise Charge, according to the manufacturer's specifications, 
can be written in the nominal working conditions as: 

ENC = a + b-C = 180e~ + 7.5e"- (7(pF) (1) 

Since a minimum ionizing particle creates on the average 27500 e~ in a 300 /im 
thick silicon layer, from eq. 1 we expect for the signal-to-noise ratio S/N on both 
sides values of this order: {S/N)x — 83; {S/N)y — 57 . Actually from 
the acquisition taken on beam test data we have obtained (fig. 5): {S/N)x — 
54 ; {S/N)y — 23 . Therefore the VAl chip in the real working condition shows 
a lower performance: this can be explained remembering that we are operating 
the preamplifier with reduced power consumption (1 mW per channel instead of 
1.3 mW per channel) by properly changing its bias conditions. We were forced to 
operate in this way because of the constraints imposed by the satellite. Nevertheless 
the signal-to-noise ratios obtained with lower power consumption are good enough 
to comply with the requirements of spatial resolution expected for the apparatus. 
We did also a check of the change of performance for the different VAl's working 
condition by measuring the noise obtained from the chip without any load attached 
(i.e. the a parameter in eq. 1 has been determined): it resulted 232 e^ instead of the 
nominal value 180 e~. 



3 Reduction and analysis of tracking system data 

36 864 electronic channels are present in the whole tracker. If an average estimated 
trigger rate of 12 Hz [11] is assumed for in-orbit acquisitions (including spuri- 
ous events, such as those related to interactions in external dead materials) a total 
amount of more than 50 GB/day can be inferred only for the spectrometer. This 
value is larger than the maximum available bandwidth for data transmission to the 
ground stations (10 GB/day). A compression procedure is then applied on-line to 
the raw data: it consists of a Zero Order Predictor-like algorithm, completed by a 
cluster finder which preserves the set of signals on adjacent channels correspond- 
ing to the passage of an ionizing particle. The amount of data to be transmitted is 
then reduced to about 5% of the initial size without degrading the spatial resolution 
of the instrument. In the off-line analysis values of non-transmitted channels are 
reconstructed on the basis of the transmitted ones. 

Data acquired on satellite consist of calibration runs, periodically repeated, and 
physics runs. In the first ones (acquisitions without particles) the pedestal (FED) 
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Fig. 3. Distribution of the R = {S/N)strip variable for all the read-out strips 
on the junction side of a ladder, as obtained from MIP tracks, orthogonal to the 
detector, on a test beam. The peak centred in zero is the noise distribution, while 
particle signals are contained in the high-5/A^ tail. The threshold Cd = 10 se- 
lects the strip whose signal is expected to be the maximum of the cluster, while 
Cind = 5 determines whether lateral strips have to be included in it. 

and the intrinsic noise (A^) of each channel are evaluated. In the physics runs the 
pedestal values are used to obtain, by a recursive procedure, the common noise 
(CN) of each VAl chip, different from an event to another. The true signal S* of a 
given channel is then obtained from the 12-bit ADC value (ADC) as: 



S = ADC - PED - CN 



When all the signals of all the electronic channels are available, those correspond- 
ing to the passage of an ionizing particle are selected through the value of the 
variable: 

R = S/N 



The R variable represents the signal-to-noise ratio of a given strip for the consid- 
ered event. A i?-value greater than a chosen threshold Cd reveals the so-called 
cluster seed; channels adjacent to the maximum on both sides are then included in 
the cluster as far as their signals fall below a second, lower threshold C^d- When a 
channel has a signal below the inclusion level C^d, search for further strips stops 
on that side of the cluster with respect to the maximum. The values Cd = 10 and 
Cind = 5 have been used for this analysis; their positions in the distribution of the 
R variable are shown in fig. 3 for the junction side of a typical detector. The number 
of strips in a cluster, according to the above inclusion rule, is called its multiplicity; 
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Fig. 4. Percentage of events with given multiplicity on both sides of the sensors 
(orthogonal MIP tracks from a test beam). 




Fig. 5. Distribution of the signal-to-noise ratio of the cluster for both X and Y 
sides of a ladder (orthogonal MIP tracks from a test beam). 



the distribution of this variable on both sides is given in fig. 4 for tracks orthogonal 
to the sensors from a test beam. Following ref. [10], the signal-to-noise ratio of the 
cluster is defined as the sum of the S/N ratios of all its strips: the corresponding 
distribution is given in fig. 5 for the same set of events as in fig. 4. 
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Fig. 6. Simulated track by a MIP in the silicon detector with emission of a 5 ray. 
The event is seen on two planes orthogonal to the strips of junction (left) and 
ohmic (right) sides. The black rectangles represent the read-out strips; the blank 
ones are the floating strips of the X, junction view. 

4 Detector simulation 

The developed simulation includes the complete chain of physical processes pro- 
duced in consequence of a crossing particle, from ionization of silicon to the dig- 
itized output signals. The first step is achieved through GEANT [16], which re- 
produces the charge generation along the track of the ionizing particle. In fig. 6 
a simulated track originating from a MIP in the silicon detector is shown for two 
orthogonal cross sections: each point along the track represents a charge packet 
and its area is proportional to the energy loss. Tracking is carried out using finite 
steps, whose maximum length has been set to be 10 fxm (corresponding to at least 
30 charge packets per track inside the 300 /im thick detector) to exploit a fine gran- 
ularity in the energy deposition. 

In a second phase the charge packets created in the sensitive volume are collected 
on the strips: their behaviour can be described as due to both drift and diffusion. 
The drift velocity v of charge carriers is given by the product of the mobility i^i and 
the electric field E: v = ■ E . Moreover the charge carriers diffuse in the detector 
material: the point-like charge packets widen according to a Gaussian distribution. 
As a consequence, the standard deviation of the packet's spatial distribution at time 
t after its generation can be written as 

a = V2Di (2) 

D being the diffusion coefficient, which depends mainly on temperature and char- 
acteristics of the material [17]. All the charge packets, located along the track, are 
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Fig. 7. Behaviour of equipotential lines (solid) and electric field lines (dashed) in 
proximity of junction {top) and ohmic {bottom) sides of the sensor. 

collected on the strips and their spatial distributions, that have different widths, con- 
tribute to the creation of a cumulative packet on the surface of the detector. To give 
an idea of the transversal size of the charge cloud, in our configuration the standard 
deviation of a packet that has covered half detector thickness is about 5 /im and it 
gets the same value for electrons and holes [18]. 

The simulation of the drift process is based on the knowledge of the electric field 
E inside the detector. Along a coordinate {z in our case) orthogonal to the p-n 
junction E can be simply obtained from the solution of the 1-dimensional Laplace's 
equation, as done in ref. [19]. In case of overdepleted detectors (i.e. if the bias 
voltage Vhias is greater than the depletion voltage Vdepi) the electric field intensity 
can be written at a distance z from the junction as: 

Vdepl I , N Vbias — Vdepl 



E{z) = -2^{d-z) 



d 



where d is the thickness of the detector, corresponding to the depletion zone. 

In the real case, the 3-dimensional map of the electric field is in principle necessary 
to determine the path in silicon of the charge carriers. Thanks to the symmetry of 
the system, this map has been obtained by solving by a numerical method the 2- 
dim. Laplace's equation for each side of the sensor [17]. In fig. 7 the equipotential 
lines and the electric field lines are shown for regions close to X and Y sides of the 
sensor, corresponding to 1/10 of the full detector thickness. For the junction side 
(top) the E-lines are orthogonal to the surface of the sensor, except in a thin layer 
(~ 10 /xm) close to the junction. Since the electric field lines are essentially straight 
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Fig. 8. Mechanism of division of a Gaussian charge packet between two subse- 
quent strips: each one collects a charge signal corresponding to the fraction of 
area closer to itself. The presence of the interstrip capacitance Cis gives rise to a 
rearrangement of the charge: a fraction of each signal appears on the output of the 
adjacent channels. As an example signals before (6) and after (a) the introduction 
of this capacitive coupling are depicted in figure. 

lines except in a small zone close to the collection plane, the drift of each charge 
packet is reproduced in the simulation simply as a "translation" of it towards the 
strips and the diffusion is considered as a Gaussian enlargement of the packet up 
to the final width given by eq. 2, where t is now the collection time of the packet. 
Finally, this charge packet is divided between adjacent strips, as illustrated in fig. 8. 

On the ohmic side (bottom) the calculation has been performed considering the 
effect of the p blocking strip. According to measurements reported in literature, 
when the junction is reverse biased such structures reach at equilibrium a potential 
which mainly depends on both the strip geometry and the doping level. Values 10- 
15 V lower than those of the adjacent ?7,+ strips have been reported [20,21]. If the p- 
stop potential is introduced as boundary condition for the solution of the Laplace's 
equation, the field lines assume a shape similar to that shown in the second graph 
of fig. 7: around the p-stop strips, a zone with a lower density of field lines can 
be observed. When particles cross the sensor in this region, charge carriers spend 
more time in the material before being collected on the n"*" strips: their diffusion 
widths increase consequently. In order to reproduce this effect in the simulation, 
the width of the resulting diffusion cloud has been parameterized for the ohmic 
side as a function of the interstrip position. 

Capacitive couplings between adjacent strips are also taken into account: a fraction 
of the signal appearing on each channel is actually visible on the adjacent output 
channels (fig. 8). In a first-order approximation, this fraction corresponds to the 
ratio r = Cis/Cdec between the first-neighbour interstrip and the decoupling 
capacitances [10]. A second-neighbour coupling has also been introduced on the 
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junction side for a fine tuning with data. Values of the fraction of transferred charge 
for both junction and ohmic sides have been obtained by tuning simulation and 
data. On the junction side an additional difficulty is the presence of the intermedi- 
ate, floating strip: charge collected by it is partly lost to the ohmic side, while the 
surviving fraction is shared upon the adjacent strips. Also this effect has been tuned 
on data. 

The electronic noise is the last physical effect to be reproduced in the simulation. To 
this purpose Gaussian fluctuations have been superimposed to the signals arising in 
the simulated electronic channels: the corresponding standard deviations have been 
extracted from a data file acquired on a test beam. 



5 Effects of a magnetic field 

The presence of a magnetic field in the region where detectors operate influences 
the collection of the charge on the strips. Assume to have a uniform magnetic field 
parallel to the strip direction and a particle beam orthogonal to the sensor. Due 
to the Lorentz force on the charge carriers, electrons and holes produced along the 
track of a crossing particle are collected at an angle (Lorentz angle) with respect 
to the track, which in the International System can be evaluated as: 

where B is the value of the magnetic induction and fifj is the so-called Hall mo- 
bility [18]. In our configuration the Lorentz angle is about 1 deg. if the value 
Hh = 0.0310 m^/Vs [18] is used along with a mean value of the magnetic field 
of 0.45 T. The vector B is orthogonal to the Y strips in PAMELA, giving rise to 
a Lorentz angle in the XZ plane (and therefore affecting the measured position 
on the junction side only). In this plane an effective angle i^^f j = tD + should 
be considered, instead of the incidence angle 'd. The magnetic field has not been 
considered in the simulation because our main present interest is the study of the 
intrinsic detector performances. Nevertheless the magnetic effect can be easily in- 
troduced in the existing code, if the characteristics of the whole spectrometer, in 
presence of straight and inclined tracks, will be investigated. 



6 Tuning of the simulation on data 

A comparison between simulation and data coming from a test beam with MIP 
tracks orthogonal to the detector (200 GeV protons gathered at CERN-SPS) en- 
ables us to check the correctness of our work. In particular the cluster charge distri- 
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Fig. 9. Cluster charge distribution (in ADC counts) on the junction side for data 
(points) and simulation (line). A good agreement is observed in both linear (left) 
and logaiithmic (right) scales. 
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Fig. 10. Comparison between the average shape of the cluster obtained from beam 
test data (points) and the corresponding simulated quantity (line) for both junction 
and ohmic sides (tracks orthogonal to the sensor). The data sample consists of 
about 20 000 clusters. 

bution, the average shape of the cluster and the r] function for data and simulation 
have been used for a fine tuning of some parameters used in the simulation. 

The distribution of the cluster charge, defined on the basis of the inclusion rules 
described in sect. 3, is reported in fig. 9. Provided charge is not lost outside the 
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Fig. 11. The r] distribution on both sides, for data (points) and simulation (line). 
In this data run tracks were approximately perpendicular to the surface of the 
detector; actually the asymmetry of the junction side distribution is due to a small 
inclination of the beam (~ 0.3 deg.), reproduced in the simulation. 



cluster, this quantity is proportional to the energy loss in the detector thickness. 
In fact, histograms in fig. 9 follow a Landau distribution, as expected for MIP's 
crossing thin layers of material. The same agreement between data and simulation, 
shown in figure for the junction side, is observed on the ohmic view. The average 
shape of the cluster is shown in fig. 10: each graph has on the abscissa the strip 
number around the maximum of the cluster and on the ordinate the fraction of sig- 
nal collected by the corresponding strip. All the clusters (about 20000 in the data 
sample) are superimposed in the histogram to obtain an average behaviour. A left- 
right asymmetry (at 3a level) can be observed in figure on the junction side for 
the strips adjacent to the maximum: this difference could be explained as due to 
a small shift (~ 1 /im) of the mask used to implant the strips in the silicon bulk 
with respect to the mask used to define the metallic strips. As a consequence, the 
capacitive coupling of a strip to the left neighbour is different from the coupling to 
the right one. This possibility has been confirmed by the manufacturer. The sim- 
ulation shown in fig. 10 includes such asymmetry. The good agreement between 
data and simulation represents a first significant test of the general model, since the 
average shape of the cluster is influenced by all the mechanisms of charge division 
and rearrangement on the strips (diffusion of charge carriers, capacitive couplings, 
signal-to-noise ratio). 



The 7] variable [10,19] is the basic parameter for our following analysis. If the pair 
of adjacent channels showing the highest signals are considered in each cluster, r] 
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can be defined as: 



where L and R are respectively the signals of the left and right channel in the pair. 
The ?7 variable can be thought as an average (weighted by the signals L and R ) 
of the positions of two strips located in and 1: if 77 ~ or 77 ~ 1, it means 
that practically only one strip collects all the cluster signal. When this variable 
is reported on graph for all the clusters generated by a set of orthogonal particles 
uniformly distributed over the sensor, the characteristic behaviours shown in fig. 1 1 
can be observed for junction and ohmic sides. This figure contains a comparison 
between the distributions of the 77 variable for simulation (line) and beam test data 
(points). Peaks near and 1 correspond to the presence of the read-out strips, while 
the peak near 0.5 on the junction side corresponds to the intermediate floating strip: 
channels adjacent to such a strip collect approximately the same amount of charge 
(i.e. L ~ i?) when the floating strip has been hit by a particle. 

Connecting the peculiar shape of the rj distribution with the detector's character- 
istics is very instructive (see ref. [10]). The positions of the lateral peaks in the 
distributions are related to the capacitive coupling of a strip to the adjacent ones: 
larger the coupling, further the peaks are from and 1, since the fraction of trans- 
ferred charge becomes more and more considerable. The peak widths are related to 
the size of the diffused packet and to the S/N ratio. Concerning the central peak 
on the X-side distribution, both the electronic noise and the loss of charge to the 
ohmic side (that resulted about 20% from a tuning with data) influence its width. 

Since the ?7-value is related to the main geometrical and electrical parameters of 
the detectors, it is a good probe to check if the code correctly reproduces the real 
detector response. A reasonable agreement between the rj distributions for data 
and simulation has been obtained by a fine-tuning of free parameters (fig. 11). 
On the basis of this positive comparison, we are confident to use the simulated 77 
distribution in the procedure of impact point reconstruction, as explained in the 
following section. 



7 Position finding algorithms 



A procedure that enables us to obtain the position of the incident particle on the 
sensor, once the signals induced over all the electronic channels are known, is called 
position finding algorithm. Such a method usually applies to clusters (identified on 
the basis of properly defined selection criteria) and allows to extract from them 
the particle impact point position. The general properties of algorithms for position 
measurements are well described in ref. [22,23]. A thorough description of the 
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Fig. 12. Simulation of two tracks inside the silicon detector and corresponding 
signals on the junction side. The typical clusters, obtained from tracks orthogo- 
nal to the detector, are shown. Channels are included in the clusters (grey bars) 
provided they satisfy the selection criteria described in sect. 3. On the left, the 
particle passes close to a read-out strip (black rectangle) which collects all the 
charge created in silicon: the corresponding cluster multiplicity is one. On the 
right, the particle crosses the detector near a floating strip (blank rectangle): the 
charge collected by it appears on the adjacent read-out channels and the corre- 
sponding cluster multiplicity is two. 

algorithms used for silicon microstrip detectors is given in ref. [10]; here we review 
only those our analysis is concerned with. 

> In the very simple digital algorithm the position of the strip which exhibits the 
biggest signal is assumed to be the incidence point of the particle. 

> Another intuitive algorithm is the Centre Of Gravity ( COG) : the particle position 
can be estimated as the weighted average of the positions Xi of the m strips 
included in the cluster: 

^ = ^^ (4) 

the weights 5*, being the corresponding signals. 

> In case of orthogonal to the sensor or slightly inclined tracks one can think of 
limiting the previous sum to the pair of adjacent strips which collect the biggest 
signals. This is a reasonable approximation because, for small incidence angles, 
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only one or two strips collect the whole charge generated in silicon by an ionizing 
particle. This is illustrated in fig. 12, where typical clusters for the junction side, 
corresponding to orthogonal tracks, are shown. Our detectors, on the X view, 
exhibit clusters whose multiplicity m is one or two for about 98% of events in 
case of perpendicular tracks; higher multiplicity is associated with emission of 
energetic 6 rays in silicon. If the events with m = 1 are separated from those with 
m = 2, the first ones are characterized by a worse accuracy in the reconstructed 
position, since in this case the information consists only in the position of the hit 
strip. 

In the assumption that two strips collect the whole particle signal, one can 
write down the following relationship (linear rj algorithm), which exploits the 
definition given in eq. 3, and corresponds to eq. 4 with the sum limited to two 
terms only: 

X = Xl + Tj ■ P 

Here xl is the position of the left strip in the pair and P is the read-out pitch. 
> Unfortunately, due to the non-uniformity of the distribution of rj, it can be proved 
that the position reconstructed by the linear 77 algorithm is systematically shifted 
with respect to the true one [18,22]. The systematically correct particle position 
is obtained by a non-linear r] algorithm containing a strictly increasing function 



The /(?7) function which produces the systematically correct impact point is the 
cumulative probability distribution function of 77; it can be estimated as: 



dN/df] being the experimental r] distribution (fig. 11), obtained by a uniform 
lightening of the sensor [10]. This method, known in literature simply as the 
7] algorithm, has been extensively employed in the PAMELA data analysis: it 
allows to reach spatial resolutions of 3 yum (junction side) and 12 fim (ohmic 
side) for orthogonal tracks [2]. In our analysis the 77 algorithm is applied if at 
least one of the channels adjacent to the maximum has a positive signal above 
the pedestal level; otherwise, the digital algorithm is used ^ . 
> In case of inclined tracks, the charge created in silicon can be collected by more 
than one or two strips, depending on the detector configuration (fig. 13). For 
some strip geometries (like that in the left picture in figure) it is natural to expect 
an improvement to the spatial resolution by an algorithm involving more than 
two strips in the determination of the impact point. For this reason, extensions 

^ Please be careful not to mistake the multiplicity of the cluster, defined in sect. 3, for the 
presence of positive channeis around its maximum, requested to appiy the afgorithm. 



(5) 




(6) 
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Fig. 13. Simulation of an inclined track inside the silicon detector and correspond- 
ing clusters for junction (left) and ohmic (right) sides. Due to the different strip 
geometries, similar incidence angles result in different cluster multiplicities. 

of the non-linear 77 algorithm can be thought, based on a larger number of strips. 
In particular, if three strips around the maximum are considered, a new variable 
can be introduced by exploiting the respective signals So, Si, S2'- 

^ S1 + 2S2 
5*0 + 5*1 + 5*2 

Similarly to what observed for r] (sect. 6), this new variable represents the weighted 
average of the positions of three strips, located in 0, 1, 2; they are chosen so that 
the central one contains the maximum signal of the cluster. A further variable, 
based on the signals collected in four channels, can be conceived in a similar 
way: 

^ gl + 2 ^2 + 3 ^3 

So -\- Si -\- S2 -\- S3 

In this case the strips are chosen so that 5*1, S2 are the biggest adjacent signals 
(those used in the rj algorithm) and 5*0, 5*3 are signals belonging to additional 
lateral strips (respectively on the left and right side of Si and S2). According 
to their definitions, 773 varies in the range [0.5; 1.5], while 774 is contained in the 
interval [1; 2]. A non-linear procedure, similar to that used for 77, can be applied 
to improve the characteristics of this position finding algorithm extended to three 
or four strips: two functions /3(ry3) and /4(r/4) can be defined on the basis of the 
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experimental distributions dN/drj^ and dN/drj^ as already done in eq. 6. The 
particle impact point will be consequently reconstructed as: 



X = Xl + f3{V3) ■ P 

x = xl + hivi) ■ P 



by the 773 algorithm 
by the 774 algorithm 



xl being the left strip in each group of three or four adjacent strips. 

Methods of impact point reconstruction for large incidence angles are also de- 
scribed in literature [10]; they are not treated here because of the limited angular 
acceptance of the PAMELA telescope, which corresponds to about [—20; 20] deg. 
for both coordinates. For such angles we expect in our strip geometry a cluster 
multiplicity of about three on the junction side and two on the ohmic side (fig. 13); 
fluctuations in these numbers, for similar incidence angles, can be produced by 
different positions of the track with respect to the strips. 

The methods presented in detail in the previous summary will be applied to the 
analysis of simulated events in the next sections. 



8 Position reconstructed by the 77 algorithm: comparison with data 

Before applying the position finding algorithms, described in the previous section, 
to the analysis of simulated events, a comparison with data concerning the position 
reconstructed by the 77 algorithm is suitable. 100 GeV electrons, with incoming 
direction orthogonal to the sensors, have been acquired without magnet on a test 
beam by a prototype telescope: it was composed of five planes, located, along the 
beam direction, in the origin, in ±145 mm and in ±165 mm. Among the gathered 
events, only those containing one cluster per view have been selected. In this case 
for each crossing particle five points {xi, Ui) can be found by the 77 algorithm, corre- 
sponding to the best estimations of the impact coordinates on every plane. A linear 
fit on these points identifies five other points (xj, yi) on the planes; the differences 
{xi — Xi) between reconstructed and interpolated positions on the junction side are 
reported for every detector in the histograms of fig. 14 (thin line). In the simula- 
tion the experimental configuration has been reproduced, including the contribution 
from the multiple scattering on every silicon plane. The same differences with re- 
spect to a fitted straight line are reported in the histograms (thick line), as done for 
data. The good agreement between data and simulation is clearly visible. 

Let us notice that the RMS's of the distributions of two different planes are not com- 
parable, since the detectors are not equidistant. Assuming the same uncertainty a 
for each reconstructed coordinate Xi, we can derive (as error propagation) the RMS 
of the distributions of {xi — Xi). By a reverse procedure, from the five RMS values 
of fig. 14 the spatial resolution of the detectors, defined as the uncertainty on the 
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Fig. 14. Distribution of the differences (xj — Xi) 
calculated for data (thin line) and simulation 
(thick line) with orthogonal tracks in absence 
of magnetic field. The experimental set-up con- 
sisted of five non-uniformly spaced planes, ar- 
ranged in a prototype telescope that acquired data 
on a test beam. 
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Fig. 15. Distributions of residuals for the rj algorithm on junction (top) and ohmic (bot- 
tom) sides. The incidence angles of the particles are, from the left column to the right 
one, 0, 10, 20 deg. 

reconstructed coordinates, can be obtained: from this data sample a = 3.0 /im [2] 
can be found. 



9 Spatial resolution for orthogonal and inclined tracks 



The aim of this section is the study of the optimal position finding algorithm for 
our detectors, as obtained from the simulation of tracks incident at various angles. 
For this purpose minimum ionizing particles have been generated with incidence 
angles uniformly distributed in the [—20; 20] deg. range. This interval has been 
subdivided in 1 deg. wide bins and consequently the simulated events have been 
classified on the basis of their incidence angles. The r] algorithm has been initially 
used on the whole angular range; a different /{t]) function (eq. 6) has been built for 
each angular bin starting from the t] distribution of the corresponding set and the 
particle position for every event has been finally reconstructed by exploiting eq. 5. 

In order to compare different algorithms, the distributions of the spatial residuals 
have been used: the residual is defined as the difference between the real particle 
impact point at half detector thickness, given by the simulation, and the recon- 
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0.5 


22.9 % 
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1.6 % 


0.4 % 


75.5 % 


23.9 % 


0.6 % 


10 ± 


0.5 


6.6 % 


84.0 % 


8.8 % 


0.6 % 


61.4 % 


37.8 % 


0.8 % 


15 ± 


0.5 


1.3 % 


65.4 % 


32.3 % 


1.0 % 


43.4 % 


55.4 % 


1.2 % 


20 ± 


0.5 


0.4 % 


30.7 % 


66.4 % 


2.5 % 


24.2 % 


72.2 % 


3.6 % 



Table 2 

Percentage of events with given multiplicity m as a function of the incidence angle. 

structed position, obtained by the chosen algorithm. For a quantitative comparison 
of different methods the RMS's of the distributions are evaluated. In fig. 15 the 
distributions of the residuals, obtained from the r\ algorithm on both junction and 
ohmic sides for 0, 10, 20 deg. incidence angles, are shown. By simply comparing 
the distribution widths, it can be noticed that at small angles the junction side ex- 
hibits a better resolution, as expected, but at 20 deg. the ohmic-side distribution is 
narrower: this is an indication of a loss of signal on the junction side out of the pair 
of strips the 77 algorithm uses, since the multiplicity of the events becomes on the 
average greater than two when the track inclination increases (in tab. 2 the percent- 
ages of events with given multiplicities m are reported for some incidence angles). 
We can expect that algorithms such as 773 or 774 will improve the resolution in this 
case. The non-Gaussian behaviour of the distributions (especially for small angles) 
can be observed too in fig. 15. On the junction side an accurate fit is provided for 
tracks perpendicular to the sensor by a generalized Lorentz distribution L (fig. 16): 



Pi ■ 



' 2 ■ {x-p2y 

P3 



Pi 



(7) 



The RMS of the distribution of residuals corresponds to the uncertainty in the re- 
constructed coordinates and consequentely can be defined as the spatial resolution 
of the detectors. The RMS of the distribution of fig. 16 (2.87 /im) can be compared 
with the spatial resolution (3.0 /im) obtained for the real detectors from the RMS's 
of the distributions reported in fig. 14, as explained in sect. 8. 

On the ohmic side, where the pitch is larger, the evident tails in fig. 15 in case of 
orthogonal tracks have been observed to be populated by some of the events in 
which the particle crosses the silicon detector some micrometers away from a strip 
and that collects all the released charge (in the Y-side geometry this may happen 
even if the distance of the particle impact point from a strip is 10 — 15 /im). In 
this case the left strip and the right one with respect to the maximum, that do not 
collect charge, can be included in the t] calculation (that always requires two strips) 
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Fig. 16. Simulated distribution of residuals on the junction side for or- 
thogonal tracks. Parameters Pi, ■ ■ ■ ,Pi are obtained from a fit with a 
generalized Lorentz distribution (eq. 7). 



with equal probability, due to noise fluctuations in these channels. If the strip away 
from the real impact point is randomly chosen between the channels adjacent to 
the maximum, large errors in the reconstructed position will be observed. It can be 
easily verified that for such events the digital algorithm is comparable to 77 [17]. 

An independent contribution to the tails of the distributions on both sides come 
from the emission of 6 rays in the detector (like that already shown in fig. 6). 
Another example of such an occurrence is shown in fig. 17 along with the cor- 
responding cluster signals: in this case the 6 ray is so energetic that two different 
peaks appear in the cluster, the highest one being produced by the Bragg's peak of 
the secondary electron. The resulting error in the reconstructed position would be 
about 150 fim for this event, if the rj algorithm were applied around the maximum. 
In spite of this large error, such a cluster can be recovered in the final data analysis 
of the spectrometer, if two possible impact points are attributed to it, each corre- 
sponding to a cluster peak: the real position can be decided by comparison with the 
impact points on the remaining detecting planes during the tracking phase. 

The behaviour of the spatial resolution as a function of the particle's incidence angle 
i!) is now studied for a sample of tracks from which pathological events (like those 
above described, with emission of energetic 5 rays) are excluded. Moreover a se- 
lection in the cluster multiplicity, depending on the track inclination, is introduced. 
The selection criteria are quite arbitrary, but they are employed only to produce 
a homogeneous sample of events where different algorithms are tested and com- 
pared. On the basis of the information contained in tab. 2, the following clusters 
have been included in our analysis: 
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Fig. 17. Track by a MIP with emission of a very energetic (about 
260 keV) 6 ray and corresponding cluster on the junction side, 
which exhibits two distinct peaks. The maximum signal corre- 
sponds to the Bragg 's peak of the secondary electron. 

Junction side (bending view) 
m = l,2 if^?<15deg. 
m = 1, 2, 3 otherwise 

Ohmic side 
m = 1,2 for every 

The 7], ?73, ?74 and COG algorithms have been applied in every bin of incidence 
angle on the junction side; on the ohmic side only i] has been utilized, because of 
the small fraction of events whose multiplicity is greater than two. 

Let us now specify an aspect of the procedure that can give rise to some misun- 
derstanding. When the t] algorithm is applied to clusters whose multiplicity is one, 
the charge released in silicon may be really collected by only one strip (the maxi- 
mum of the cluster) and both the signals of the adjacent channels may be negative 
for noise fluctuations. In this case only the signal of the maximum is used in the 
impact point reconstruction and r] reduces to the digital algorithm. No restrictions 
are therefore introduced to such clusters, whose information is fully contained in 
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Fig. 18. Spatial resolution obtained by some position finding algorithms as a function of 
the incidence angle on both junction and ohmic sides. 

the position of the maximum. The same criteria are followed for the other methods: 
as a general rule, when a sufficient number of strips with positive signals are not 
present in a cluster to apply a given algorithm, a procedure involving less channels 
is used instead, according to the sequence: 

V4: ^ Vs ^ V ^ digital algorithm. 

In fig. 18 the results provided by the mentioned methods on both sides are shown. 
On the junction side it can be noticed that the r] algorithm gives the best resolution 
for orthogonal or slightly inclined tracks, as already stated in literature [10,18,19]: 
in our configuration t] is the preferable method up to about 10 deg. For larger angles, 
the original methods 773 and 774 give a more accurate particle position, compared to 
the Centre Of Gravity, suggested in the ref. [10] for intermediate incidence angles. 
In particular, 773 can be used from about 10 to 15 deg. and 774 from 15 to 20 deg. The 
spatial resolution ranges between less than 3 //m and about 7.5 /im on the angular 
interval corresponding to the acceptance of the spectrometer. On the ohmic side a 
better resolution (about 8 /im RMS) is achieved for inclined (t? ~ 10 deg.) tracks. 
This is due to a greater percentage of clusters with multiplicity m = 2 compared to 
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case of perpendicular tracks, which result in a higher accuracy of the reconstructed 
position. 



10 The next step of this work 

Let us now briefly discuss how the spatial resolution, achieved at a given angle, 
influences the scientific objectives of the PAMELA experiment. The silicon mi- 
crostrip detectors are the sensitive elements of a magnetic spectrometer that deter- 
mines the momentum of the crossing charged particles. The obtained uncertainty 
Ap/p in the reconstructed momentum is directly related, besides the magnetic 
field value and the geometry of the spectrometer, to the spatial resolution on the 
X (bending) view. Usually the performances of the instrument are given in terms 
of the rigidity r, defined as momentum-to-charge ratio: when Ar/r = 100 % the 
spectrometer has reached the Maximum Detectable Rigidity (MDR). In R^.MELA 
the value MDR = (1183 ± 54) GV/c [24] has been measured on a test beam for 
electrons whose arrival directions are orthogonal to the sensors. This value is in- 
fluenced by the presence of tails in the distribution of residuals that can in princi- 
ple be associated with errors in the charge sign of the reconstructed events. This 
problem becomes relevant in R^MELA because few antiparticles have to be iden- 
tified in a large particle background (the antiparticle/particle ratio is of the order 
of 10^"^). This drawback is known as spillover and actually limits the measure- 
ments of antiparticle abundances in R^MELA to about 190 GeV/c for antiprotons 
and 270 GeV/c for positrons [2], in spite of a much larger MDR. For these reasons 
the data analysis needs accurate procedures that do not disperse the peculiar char- 
acteristics of the detectors. Results from this simulation will be therefore exploited 
to obtain the best algorithm of position reconstruction depending on the particle 
incidence angle and on the cluster multiplicity. Moreover, the knowledge of the 
width of the distribution of residuals in different conditions enables us to "weight" 
each impact point determination while tracking a particle through the spectrometer. 
For instance, a high-multiplicity cluster, occurring in a quasi-orthogonal track, is 
probably connected with a 5 ray emission: so the corresponding position should 
reasonably have a lower weight in the track fit. If such a correlation is considered, 
the accuracy in the momentum reconstruction will be significantly improved. 



11 Conclusions 

The simulation described in this paper takes into account the characteristics of 
the double-sided silicon microstrip detectors of the R^.MELA spectrometer and 
correctly reproduces their observed performances. In particular the simulation has 
been applied to the analysis of the best obtainable spatial resolution as a function 
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of the incidence angle of the incoming particle. New position finding algorithms 
(?73 and 774) have been successfully introduced to improve the resolution for non- 
orthogonal tracks on the bending view. The obtained distributions of the spatial 
residuals at different angles and multiplicities can be used in a simulation of the 
whole tracking system to study the particle contamination in antiparticle measure- 
ments and to calibrate in the best way the tracking procedure. 
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